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^'- 

'~~^ ' Direct and indirect dark matter searches can potentially lead to the discovery of Su- 

C^^ . persymmetry. In the specific context of heterotic orbifold conipactification scenarios, it is 
even possible to get some clues on the SUSY breaking mechanism of the more fundamen- 

s^ . tal underlying theory. In this paper, we investigate the prospects for indirect detection of 

*yv ' neutralino dark matter via gamma-ray (continuum and line) and synchrotron radiation from 

QQ . the Galactic center, in the context of AMSB scenarios in an effective heterotic framework, 

^^ ' which constitutes a consistent SUSY model between high energy theory and low energy 

>0 , phenomenology. 

^: 

^ ■ I. INTRODUCTION 

I I Many independent measurements provide convincing evidence that most of the matter in the 

flj ' Universe is dark, i.e. non-baryonic and of unknown nature, the most common dark matter (DM) 

^ ■ candidate being a Weakly-Interacting Massive Particle (WIMP) jj[|. Supersymmetric (SUSY) 

^ i extensions of the Standard Model (SM) naturally predict a massive weakly interacting particle 

(WIMP) called neutralino {xi = x)j which, in most versions of SUSY models, is the lightest 

supersymmetric particle (LSP), and is made stable by virtue of the conservation of R-parity. 

Dark matter searches are often discussed in the very popular mSUGRA, a.k.a CMSSM (con- 
strained MSSM), model ^j where one assumes unification of the soft SUSY breaking parameters 
at the GUT scale. This assumption is essentially a reduction of the parameter space and has no 
real theoretical motivation. More refined SUSY models have been studied in the context of DM 
searches, relaxing the hypothesis of universality 0, 0, H 0, 0, S, II ED, E3, 113, 0, Q, IH, El • 



X 



Here we focus on effective string inspired models, which represent a possible consistent frame- 
work linking fundamental high energy theory with low energy physics. Recently, the full one loop 
soft supersvmmetry breaking terms, in a large class of superstring effective theories, have been 
calculated [l3] , based on orbifold compactifications of the weakly-coupled heterotic string (includ- 
ing the so-called anomaly mediated contributions). The parameter space in this class of models 
^ las already been severely constrained taking into account accelerator and relic density constraints 



1^1, direct or indirect detection of DM from the Sun jl^, or benchmark models at the Tevatron 
20|. In this model, supersymmetry breaking is transmitted by the auxiliary ('i*") fields of the 
compactification moduli T", whose expectation values determine the size of the compact manifold, 
or by the dilaton auxiliary field S, whose vacuum expectation value determines the magnitude of 
the coupling constant ^STR- The spectrum and couplings of SUSY particles depend crucially on 
the contributions of each hidden fields to the breaking mechanism. 



There exist many ways of probing the existence of neutrahno dark matter and super symmetry. 
Of course, accelerator searches are among the most promising ones. Recent results give a mass limit 
of 45 GeV for the lightest neutralino, in the framework of mSUGRA |25|. Furthermore, interesting 
constraints have been derived from many experiments on direct and indirect detection of dark 
matter, whose sensitivities are expected to increase by orders of magnitude in the near future (see 
Ref. jll] and references therein) . Direct detection is made via the measurement of the recoil energy 
of the nuclei of the detector after an elastic scattering with a WIMP. Indirect detection experiments 
look for the products of annihilation of the LSP in the Sun, Galactic halo or external galaxies (see 
Ref. ^1] and references therein). 

It has been emphasized in [1^ that one-ton detectors are needed for direct searches to test all 
the parameter space in scenarios dominated by the dilaton F— breaking terms, and that moduli 
dominated scenarios are generically not detectable by dark matter searches (direct or indirect 
detection with neutrinos coming from the Sun) if we combine favored relic density and accelerator 
constraints. This means that such DM searches could be a way to distinguish the origin(s) of the 
SUSY breaking mechanism, and to restrict or even determine some fundamental parameters of the 
model (of stringy nature or not). 

The purpose of the present paper is to investigate DM indirect detection with high energy 
gamma rays and synchrotron radiation from neutralino annihilation at the Galactic center, in the 
class of weakly coupled heterotic string models discussed above 17i.ll8l. Il9l. 1201]. 



The paper is organized as follows: after a brief survey on the motivations, phenomenology and 
construction of the effective string models studied here, we analyze in section 2 and 3 the influence 
of the different coefficients parameterizing the SUSY breaking, on the gamma-ray and synchrotron 
emission from the Galactic center. In section 4, we analyze the prospects of observation of gamma- 
ray fluxes with present and foreseen experiments and discuss synchrotron radiation in section 5. 
We present our conclusions in section 6. 

II. THEORETICAL FRAMEWORK 
A. Structure of heterotic orbifolds models at one loop 

The task of string phenomenology is to make contact between the high energy string theory, and 
the low energy world. For this purpose, we need to build a superstring theory in four dimensions, 
able to reproduce the Standard Model gauge group, three generations of squarks, and a coherent 
mechanism of SUSY breaking. We set our analysis in the framework of orbifold compactifications 
of the heterotic string, within the context of supergravity effective theory, we focus on models where 
the action is dominated by 1-loop contributions to soft breaking terms. The key property of the 
models is the non-universality of soft terms, consequences of the beta-function appearing in the 
superconformal anomalies. This non-universality gives a peculiar phenomenology in the gaugino 
and the scalar sector, modifying considerably the predictions of mSUGRA models. In fact, these 
string-motivated models show a new behavior, that interpolates between the phenomenology of 
unified supergravity models (mSUGRA) and models dominated by the superconformal anomalies 
(AMSB). The constraints arising from accelerator searches, and some dark matter aspects (direct 
and neutrino indirect detections) have been already studied in J18l.ll9|. We extend here the analysis 
to gamma-ray and synchrotron emission from the Galactic center. 

We provide a phenomenological study within the context of orbifold compactifications of the 
weakly-coupled heterotic string, where we distinguish two regimes: 

• In the first one, the SUSY breaking is transmitted by the compactification moduli T", whose 
vacuum expectation values determine the size of the compact manifold. Generic (0,2) orbifold 



models contain three Tq, moduli fields. We considered a situation in which only an "overall 
modulus T" field contributes to the SUSY-breaking. The use of an overall modulus T is 
equivalent to the assumption that the three T^ fields of generic orbifold models give similar 
contributions to SUSY-breaking. This is expected in the absence of some dynamical effect 
that would strongly discriminate the three moduli. 

In the second case, it is the dilaton field S present in any four-dimensional string (whose 
vacuum expectation value determine the magnitude of the unified coupling constant gsTR at 
the string scale), that transmits, via its auxiliary component, the SUSY breaking. We work 
in the context of models in which string nonperturbative corrections to the Kahler potential 
act to stabilize the dilaton in the presence of gaugino condensation 2l|, |2j| . 



The origin of the soft breaking terms are completely different in the two scenarios. Some 
are coming from the superconformal anomalies and are non-universal (proportional to the beta- 
function of the Standard Model gauge groups), others are generated in the hidden sector (from 
Green-Schwarz mechanism or gaugino condensation) and are thus universal. This mixture between 
universality and non-universality gives the richness of the phenomenology in this type of effective 
string models and confirms the interest of non-universal studies in the prospect of super symmetric 
dark matter detection, non-universality being in this case connected with the basic properties of 
the model. 

B. SUSY parameters and constraints 

In all the models considered here, the lightest neutralino is the LSP. Its phenomenology (mainly 
the couplings and the spectrum) is determined by its mass matrix, 
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written in the {B,W'^,H^, H^) basis, where, B, W"^, H^, H^ represents respectively the B-ino, 
W-ino and down, up-Higgsinos fields. Mi, M2 and fi are the bino, wino, and Higgs-Higgsino mass 
parameters respectively. Tan/3 is the ratio of the vev of the two Higgs doublet fields. This matrix 
can be diagonalized by a single unitary matrix z such that we can express the LSP as 

X = Z^lB + Z^2W + Z^3Hi + Z^4^H2. (2) 

The parameter fi is obtained under the requirement that the electroweak symmetry breaking 
(EWSB) takes place. This is done by computing the complete one-loop corrected effective potential. 
The effective /i-term is calculated from the EWSB condition after minimization of the scalar 
potential 
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H, 



+ ^^l, ) - ("^k + -^"^k ) tan2 /? 1 



u = tt-t; M7 (3) 

^ tan2 p-1 2 ^ ^ ^ 

The sign of /x is not determined and left as free parameter. 

We clearly see that the nature of the neutralino (and its couplings) depends crucially on the 
hierarchy between the parameters Mi, M2 and /U. 



The moduli dominated scenario 



In the moduU dominated scenario, the one loop order supersymmetric SUSY breaking terms at 
GUT scale can be written [HSlii]: 
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(4) 
(5) 
(6) 



where Ma and Mj are the soft masses for the gauginos and scalars and Ai, the trilinear coupling. 
ba is the beta-function coefficient for the gauge group Ga'. 



ba 



167r2 



3Ca — 2_^ Ca 



(7) 



where Ca, C* are the quadratic Casimir operators for the group Ga in the adjoint representation 
and in the representation of the field i respectively. F^ and F^ are the auxiliary fields for the 
dilaton and the Kahler modulus, respectively, M is the supergravity auxiliary fields whose vacuum 
expectation value {vev) determines the gravitino mass m^/2 = ~|-^i ^^^ Sqs is the Green-Schwarz 
coefficient which is a (negative) integer between and —90. The function G2(T,T) is proportional 
to the Eisenstein function and vanishes when T is stabilized at one of its two self-dual points. 
From Eq.Q, it follows that when the moduli are stabilized at a self dual point, only the second 
term contributes to gaugino masses. This is precisely the "anomaly mediated" contribution. The 
loop contributions have been computed using the Pauli-Villars (PV) regularization procedure. 
The PV regular fields mimic the heavy string modes that regulate the full string amplitude. The 
phenomenological parameter p which represents the effective modular weight of the PV fields is 
constrained to be not larger than 1, though it can be negative in value. Thus the scalar squared 
mass for all matter fields is in general non-zero and positive at one loop (only the Higgs can have 
a negative running squared mass). The limiting case of p = 1, where the scalar masses are zero at 
one loop level and for which we recover a sequestered sector limit, occurs when the regulating PV 
fields and the mass-generating PV fields have the same dependance on the Kahler moduli. Another 
reasonable possibility is that the PV masses are independent of the moduli, in which case we would 
have p = 0. and 7^ is related to the anomalous dimension through -ff = jiSj. (see 121,113)1201 for 
notations and conventions) 

We clearly see in these formulae the competition between universal terms and non-universal 
ones. The scalar mass terms are all non-universal and proportional to their anomalous dimension 
and thus loop suppressed. The Green-Schwarz mechanism generates universal breaking terms for 
the gauginos (proportional to Sqs) whereas superconformal anomalies introduce non-universal 
contributions (proportional to ba)- The nature of the neutralino thus depends mainly on the value 
of the Green-Schwarz counterterm 6qs, whereas the mass scale is the gravitino mass m^/2. We 
illustrate in Figs ^ the neutralino mass and its relic density as a function of SqSj vertical directions 
correspond to different values of m^/2 ■ 
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FIG. 1: a) Neutralino mass b) Relic density as a function of 6gs for t — 0.25, p — 0, tan/5 = 5 in blue 
crosses (resp. 35 in red boxes) and Oi?7i3/2ilO (resp. 20) TeV. The hole for small \Sgs\ values corresponds 
to a gluino LSP corridor and the left boundary of the clouds corresponds to the stau LSP region (see Figs^ 
and lT^ 



D. The dilaton dominated scenario 



We turn now to a scenario where the dilaton is the primary source of supersymmetry breaking 
in the observable sector. It is well known that if we use the standard Kahler potential derived from 
the tree level string theory, it is very difficult to stabilize the dilaton at acceptable weak-coupling 
values. We postulate, in our study, nonperturbative correction of stringy origin to the dilaton 
Kahler potential. In that case, one condensate can stabilize the dilaton at weak coupling, while 
simultaneously ensuring vanishing expectation values at the minimum of the potential. The key 
feature of such models is the deviation of the dilaton Kahler metric from its tree level value. If 
we imagine the superpotential for the dilaton having the form W{S) oc e~^^'^+, with 6+ being 
the largest beta-function coefficient among the condensing gauge groups of the hidden sector, 
then we are led to consider the phenornenology of models given by the following pattern of soft 
supersymmetry breaking terms 
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where //uv is an ultraviolet regularization scale (of the order of the string scale Mstr) and ^u/j the 
renormalization scale (taken at the boundary value of MGUT)|5a]. Moreover 
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FIG. 2: a) Neutralino mass b) Relic density as a function of 6+ for t = 0.25, p = 0, tan/3 = 5 in blue crosses 
(resp. 35 in red boxes) and 0irn3/2il5 (resp. 4) TeV. 
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(see |l7l . Ila . |20| for notations and conventions) to ensure a vanishing vacuum energy in the dilaton- 
dominated hmit. 

The phenomenology of a dilaton-dominated scenario is completely different from a moduli- 
dominated one. If we look at formulae Q and (jlOp . it is clear that we are in a domain of heavy 
squarks and sleptons (of the order of magnitude of the gravitino mass) and relatively light gauginos. 
Indeed, the factor 6-|-, as it contains a loop factor, can suppress the magnitude of the auxiliary field 
F^ relative to that of the supergravity auxiliary field M through the relation (|11|) . The resulting 
gaugino soft breaking terms are less universal for low values of 6+. We illustrate neutralino mass 
and relic density as a function of 6+ in Figs El 



III. SUPERSYMMETRIC DARK MATTER PHENOMENOLOGY 



A. Relic density 



The relic density of neutralinos depends on their composition. In a large parameter space of 
the mSUGRA model, the bino-like nature of the lightest neutralino x implies a rather low rate of 
annihilation. 

Different processes lead to interesting neutralino relic density. Let us begin with a brief survey. 
For a bino neutralino one needs sfermion coannihilation or annihilation into the pseudo-scalar A to 



have a cosmologically favoured abundance of neutralino. If the hghtest neutrahno has a dominant 
wino component the rehc density drops because of efficient annihilations into gauge bosons as well 
as strong xxt coannihilations. For a non negligible higgsino component, the neutralino annihilates 
into gauge bosons or ti and relic density is also decreased by xxt ^^'^ XX2 coannihilations. 

In Fign(resp. HJ, decreasing 5gs (resp. b^) in the moduli (resp. dilaton) dominated scenario 
corresponds to going from a dominant wino (resp. bino) to a dominant bino (resp. higgsino) 
component. The theoretically predicted relic density can be checked against cosmological obser- 
vations. In particular, recent data on the cosmological microwave background (CMB) from the 
WMAP satellite |2a] constrain the dark matter relic density to be (at the 2a level) in the range 
ficDMh^ = 0.1126l°;oJ8J. We indicate on FigsUJa) andUJa) the WMAP favoured range, keeping in 
mind that the relic density calculation is very sensitive to SUSY mass and coupling uncertainties. 
In the following, we plot models in a generous range of relic density 0.03 < O/i^ < 0.3, indicating 
also the WMAP constraint. 



B. Dark matter profiles 

A crucial ingredient for the calculation of annihilation fluxes is the density profile of dark matter, 
which is usually parameterized as 

^^""^ " {r/Rp[l + (r/i?)"](/5-7)/« ■ (^^) 

where r is the galacto-centric coordinate, i? is a characteristic length and a, (3 and 7 are free 
parameters. Unfortunately, large uncertainties are associated with such profiles, especially in the 
innermost regions of galaxies, i.e. regions where, in many cases, most of the signal comes from. 

N-body simulations suggest the existence of 'cuspy' profiles, following a power law p{r) = r~^ 
where 7 should be ~ 1 at small radii, although its exact value is under debate. Several groups 
tried to reproduce the initial results of Navarro, Frenk & White [23|, who found 7 = 1, but reached 
different conclusions. In Tab. ^ we give the values of the parameters (a, /?, 7) for some of the most 
widely used profile models, namely the Kravtsov et al. (Kra, J28| | ), Navarro, Frenk and White 
(NFW, [13]), Moore et al. (Moore, Hi]) and modified isothermal (Iso, e.g. [13]) profiles. The 
most recent N-body simulations (Hayashi et al 2003) suggest that profiles do not approach power 
laws with a well defined index at very small radii. Profiles continue to become shallower, i.e. the 
(negative) logarithmic slope becomes higher, when moving towards the centre. 

Furthermore, the presence of a 3.6 x 10^ solar masses black hole lying at the Galactic Center (see 
e.g.Ref. I^JI) could possibly modify the profile of dark matter, that would accrete on it producing 
a so-called 'spike' J3j], leading to an enhancement of the annihilation flux by several orders of 
magnitude. The prospects of indirect detection of dark matter in presence of such a spike have 
been discussed in Ref. J3^ (see ibid, for a discussion of the dynamical effects, recently proposed 
in literature, that could potentially destroy the spike). The observational situation is even more 
unclear. Density profiles are usually reconstructed from the observation of rotation curves of 
galaxies, in particular of low surface brightness galaxies (LSB), that are thought to be dark matter 
dominated, de Blok et al. |2J] used this method to claim the inconsistency of the observed 'flat' 
profiles, with the cuspy profiles predicted by n-body simulations. Other groups 35|, |3g] claimed 
instead that cuspy profiles are compatible with observations. Hayashi et al |22| compared the 
observational data with their numerical simulations (not with fitting formulae of their simulations) 
and found no significant discrepancy in most cases. They attributed the remaining discrepancies 
to the difference between circular velocities and gas rotation speed in realistic triaxial halos. For 
more information on DM profiles see Ref. [jj. We assume in the following a NFW profile, and an 





a /3 7 


R(kpc) J{10~^) 


Kra 


2.0 3.0 0.4 


10.0 2.166 X 10^ 


NFW 


1.0 3.0 1.0 


20 1.352 X 10^ 


Moore 


1.5 3.0 1.5 


28.0 1.544 X 10^ 


Iso 


2.0 2.0 


3.5 2.868 X 10^ 



TABLE I: Parameters of some widely used density profiles models and corresponding value of J(10 '') 
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sd, the results for other profiles can be easily obtained using the numbers 



C. Flux of secondary particles 



Indirect detection of dark matter is based on the observation of secondary particles originating 
from dark matter annihilations in a cosmic storage area like the galactic halo. The study presented 
here is conceptually similar to the one in Ref. [Sg], devoted to indirect detection of Kaluza-Klein 
dark matter. 

The observed flux of secondary particle of species i, from the annihilation of dark matter particles 
of mass M and annihilation cross section av, from a direction ip and at energy E, can be expressed 
as (e.g. Ref. ^,^]) 



^i{ip,E) =av 
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where dNi/dE is the spectrum of secondary particles per annihilation and r^ 
with Rq ~ 8.5kpc is the solar distance to the galactic center. 

It is customary (see [Sfl,]), in order to separate the factors depending on astrophysics from those 
depending only on particle physics, to introduce the quantity J{'4j) 



JW 



8.5 kpc V 0.3 GeV/cm'^ / Jline of sight 



dsp (r(s,'i/')) • 



(15) 



We then define J{AQ) as the average of J(V') over a spherical region of solid angle AQ, centered on 
ip = 0. The values of J(AO = 10~^) are shown in the last column of Tab. ^ for the corresponding 
density profiles. In what follows, we assume a NFW profile : results corresponding to the other 
profiles of Table 1) can be deduced by a rescaling of J 

WIMP annihilation in the Galactic Halo produce a flux of gamma-rays, with either a continuum 
(coming from vr" decay after the shower from hadronization of the neutralino annihilation prod- 
ucts (quarks or gauge/higgs bosons)) or a monochromatic distribution (produced from one-loop 
annihilation into 77 or Zj). Monochromatic processes are loop suppressed, and contribute little 
to the 7 spectrum. But, as the neutralino can be considered at rest, the monochromatic ray is at 
an energy equal to the neutralino mass. This could be one of the most promising signals for the 
discovery of supersymmetry in dark matter searches. 

At the same time, annihilating neutralinos can also generate positron fluxes. Cosmic-ray elec- 
trons and positrons interact with the interstellar medium through synchrotron radiation and inverse 
Compton scattering. The observed e flux is dominated by primary electrons from acceleration 



sites, however, 10 percent of the total flux is made of secondary electrons and positrons produced 
(in equal number) by the interaction of the primaries with the interstellar medium. 

Data from HEAT (High Energ y A ntimatter Telescope) suggest the existence of a "bump" in 
the positron flux around 10 GeV [4i|. Halos of neutralinos could be a new source of positron 
explaining this excess, as there is no standard mechanism that would produce a signal at such a 
high energy. However, it is difficult to reproduce the normalization of the observed positron fraction 
even assuming that fluxes are "boosted" , e.g. because of the presence of dense DM substructures 
in the solar neighborhood (see the discussion and references in Ref. [l|). Moreover, these secondary 
e could propagate in the galactic magnetic field, generating synchrotron radiation (see below). 

D. Experiments 

Up to now, detectors have explored energy ranges above ~ 300 GeV (ground-based telescopes 
like Whipple or Cangaroo) or below ~ 20 GeV (EGRET). But the region which might turn to 
be the most promising one for neutralino physics is the intermediate region (30 GeV < E^ < 
300 GeV) that will be explored by the next generation of detectors. In this paper we focus on 
the complementarity of two type of telescope : a ground-based Atmospheric Gerenkov Telescope 
{High Energy Stereoscopic System, HESS) and a satellite experiment {Gamma-Ray Large Area 
Space Telescope, GLAST). 

Ground-based and satellite experiments are complementary sources of data for super symmetric 
dark matter searches. While satellite experiments allow for a lower energy threshold (less than 1 
GeV for GLAST), a better energy resolution and a longer exposure time, the small effective area 
of the order of the square meter limits the sensitivity to high energy photons. On the contrary, 
the large efi^ective areas (~ 0.1 km^ for HESS) of ground-based Gerenkov telescopes permit the 
measurement of very high energy fluxes with a higher energy threshold {E^ > 60 GeV for HESS) 

IV. ANALYSIS 

Gamma ray fiuxes strongly depend on the phenomenology of the neutralino. On the other 
hand, the nature of the neutralino is determined by the fundamental parameters of the models at 
high scale. We first analyze the dependence of neutralino annihilation on the main parameters of 
moduli dominated scenario ("1,3/2, ^gs, tan/3) or dilaton dominated scenario {m^/2^ &+) ^^^ then 
the resulting gamma fiux dependence on these parameters. We then compare the predictions of our 
models with the sensitivities of present and future gamma ray experiments. All calculations are 
achieved assuming an intermediate NEW profile: the results for other distributions can be easily 
deduced from a rescaling by using the J values in the table ^ 

A. Neutralino annihilation 

The main neutralino annihilation processes are XX ~^ tt {a^j oc [2:^3(4)]^), XX ~^ W~^W~ 
{(Tww o<^ [-2x3(4)^12]^ and/or [zx2^ii]^ where ^11(2) is the wino (higgsino) component of the ex- 
changed chargino) and xX ~^ bb {a^^j^ oc [2^i(2)'2;^3(4)]^)- Thus, both neutralino wino and higgsino 
components strongly enhance annihilation. The processes are shown in Eigs. |21 The opening of the 
different channels obviously depends on kinematics {m^ versus mb,mw,mt)- For a dominant bino 
LSP (Ml << M2, /x), the s-channel A exchange (Fig. |3l) is the only one present. This process is 
dominating for low rriA especially when m^ approaches txiaI'^ (A-pole) and is thus favored for high 
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a) b) c) 

FIG. 3: The dominant Feynman graph contributing for the process X X ^ 7 7- 



values of tan/3 thanks to the 66 couphng to the pseudoscalar. When kinematically ahowed, anni- 
hilation into it final state by Z exchange can be also open despite the coupling suppression. For a 
wino-like neutralino (M2 << Mi, ^u), the annihilation into WW final state through Xi^ exchange 
(Fig. |3t) is dominant thanks to the xxt^ coupling and its propagator factor (m^ ~ M2 ~ it^x^^- 
For an higgsino-like neutralino (/i << Mi, M2), when my, > mt, the Z exchange {ti final state) 
is enhanced through the XX^ coupling. But for even lower values of // (ie neutralino mass), ti 
final state is kinematically closed, and the neutralinos self annihilate into W^W~ or 66 (when 
rrix < rny/) final states. 

The differential energy spectrum of the photon depends a lot on the primary product of the 
neutralino annihilation. The gamma-ray spectrum was simulated in j5j] and fitted well with the 
function dN^/dx = ae-''^■/x^■^ where x = E^/rUx and (a, 6) = (0.73,7.76) for WW and ZZ, 
(1.0, 10.7) for 66, (1.1, 15.1) for ti, and (0.95,6.5) for uu. As we will see later on, final states with 
gauge boson will produce the harder spectrum. 



B. Gamma-ray flux in the moduli dominated scenario 



The integrated gamma flux gives a good understanding of the phenomenological aspects of the 
moduli sector of the theory. 

As an illustration, dominant annihilation branching ratios are represented as a function of 
5gs and 7713/2 in Figs. E] for tan/3 = 35 (and in appendix Figs. El for tan/? = 5) without any 
experimental cut on the parameter space. We also show in Fig. |2K the integrated (from IGeV to 
m^) gamma flux for a NFW halo profile as a function of the Green-Schwarz counterterm (5gs for 
two values of tan/3 (5,35) and 7713/2 = 10 TeV. To complete our illustration, we have plotted in Fig. 
Eb the mass spectrum 277i^ and Ma for the same set of parameters. The mass of the neutralino 
(increasing with |(^g5|) E<1- B]) and its nature (from wino to bino with increasing |5gs|) determine 
in a large part the integrated flux. Indeed the LSP is mainly bino-like for — 20 < 5gs < — 3 and 
wino-like|57l| for — 3 < 5gs < 0. 

Fig. inK illustrates first the general increase of the flux as a function of 5gs, due to the decrease 
of the neutralino mass (Fig. [SJa). For 5gs running from -10 to -4, we observe a higher flux for 
higher tan/? because the A-pole channel is kinematically open for 5gs ~ ~5 at tan/3 = 35 (value 
of 5gs that gives 2777^ > Ma)- This process produces a large amount of 66 final states (Fig. ^ 
thanks to the coupling Abb, proportional to tan/3. The low tan/3 (Fig. I12|) dominant channel 



is the tt, proportional to 



mtm^ 



For — 2 < 6gs < 0, the neutralino is entirely wino, leading to 



annihilation into WW (Fig. \^. Indeed, a crude approximation tell us that the WW cross section 
is proportional to -2^2 (~ 1 ^01 5gs = ~1) and the annihilation into ti, proportional to the square 
of the XX^ coupling ( [2^3 + 2^4] ~ 10^^ for Sgs = —4.). This explains the 5 orders of magnitude 
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Moduli : Branching ratios : tb=35 




W"W 



W"H' 




FIG. 4: Moduli domination regime: dominant annihilation branching ratios in the {SgSt 1^3/2) plane for 
tan/3 = 35, t — 0.25, p — 0. Regions with gluino LSP are indicated. We also show in grey the kinematically 
forbidden region for each channel. 



observed on the Fig. [5] between the small and large values of {dcsl for tan/3 = 35. 

For 6gs ~ and small values of 171^/2, when the WW channel is kinematically closed, dominant 
7-ray contribution comes from the annihilation into bb through Z and A exchange (suppressed 
because proportional to irib)- For tan/3 = 5 and for higher m^/2 values one can have ti or gg 
annihilation channels which explains the change of slope of the flux with tan /3 = 5 on Fig. |S1 



C. Gamma- ray flux in the dilaton dominated scenario 



In scenarios where the supersymmetry is spontaneously broken through the vev of the dilaton 
auxiliary field, the (3 function 6_|_ of the first gaugino condensing group plays a role similar to the 
Green-Schwarz counterterm \6gs\ in the moduli domination scenario. Since 6+ does not depend 
on the gauge group indices of the Standard Model 5*^^(3) x SU{2) x C/(l), it contributes universally 
to the gaugino mass breaking terms. Thus, the phenomenology of dilaton scenarios at large 6+ in 
the gaugino sector is similar to the one of the moduli dominated models at large | Sqs I • The scalar 
sector is responsible of most of the phenomenological differences. Indeed, the scalar masses are non- 
universal and suppressed in moduli like scenario because of the loop anomaly factors suppression 
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FIG. 5: Moduli-dominated scenario, a) The integrated gamma flux {E~f > 1 GeV) for a NFW profile in the moduli 
parameter space as a function of the Green Schwarz counterterm Sqs, for 77x3/2 = 10 TeV, t = 0.25, p ~ 0, tan /3 = 5 
(blue crosses) and 35 (red boxes) and b) the relevant mass spectrum of as a function of Sgs for the same values of 
tan p. See the text for the comments. 



7i ~ Yg^ (Eq. ini) , whereas they are universal and large in a dilaton dominated scenario (mo ~ 1^3/2 
Eq. llUj) . The main consequences are heavier pseudoscalar[5g| A and M^ breaking term at GUT 
scale, that drives the /i parameter to lower values, especially in high tan/3 regimes, Eq.(I2J. 

As in the moduli dominated scenario, we show the dominant annihilation branching ratios, 
in the (6+, 7713/2,) plane for tan/3 = 35 in FigsEl (and in appendix in Figs^l for tan/3 = 5) to 
highlight the different annihilation channel dependence on fundamental parameter. We also show 
(in Fig. [7^) the integrated gamma flux for a threshold E^ > 1 GeV, from the Galactic Center for 
a NFW profile, as a function of b-^., for two values of tan/3 (5 and 35) and the higgsino fraction of 
the lightest neutralino in Fig. [7|3. 

The effect of 6+ on the neutralino nature is complex and indirect. It acts on M3 through the F 
contribution with a different sign contribution than the superconformal contribution proportional 
to 63 (Eq. (jHI). Thus, lowering 6_|_ acts through the M3 dependence of the RGE on M^ . Lower 6_|_ 
means lower M3, Eq.® leading to higher M^ through the RGE implying a lower fj, parameter, 
Eq.(I2I): neutralino is mainly higgsino for 6+ < 0.1(0.3) and tan/3=5 (35) as we can see in Fig. |2b- 
Though quite low, its wino component is not completely negligible. 

For a fixed value of m^/2, decreasing 6_|_ lead to a higher higgsino fraction favoring WW, ZZ 
channels and the tt one when kinematically open (i.e for high 7713/2 values) compared to the bb 
channel. For low 7773/2 values the wino component also favors the WW channel. When 777;;^ < mw, 
annihilation into bb is completely dominant, the other processes being kinematically closed. This 
is illustrated on Figs. IHlfor tan/3 = 35 and in appendix Fig. ^Jfor tan/3 = 5. For a given 7773/2, 
higher values of tan/3 give better fluxes because it enhances the higgsino fraction of the neutralino 
through Eq.(|3J. The neutralino mass increasing with 6+ (Fig|5^) explains the general evolution 
of the gamma flux (oc -\ ~ ■^) (Fig. [7^). The gamma flux follows also the higgsino fraction 

(Fig. [TIa) which enhances the annihilation. The different peaks correspond to changes in hierarchy 
between Mi, ^,Mw in the neutralino mass matrix Eq. ^ exchanging the annihilation processes 
(Fig. inj. The different flux slopes (FigCJ correspond to the different spectrum shape of each 
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FIG. 6: Dilaton domination regime: dominant annihilation branching ratios in the (6_|_, 7713/2) plane for 
tan/3 = 35, t = 0.25, p — 0. Regions where radiative electroweak symmetry breaking can not occur are 
indicated (No REWSB). We also show in grey the kinematically forbidden region for each channel. 



neutralino annihilation final states (FigEl). 



D. Comparison of models in view of the experimental sensitivities 



The phenomenology of the two SUSY breaking mediation sectors strongly depends on the 
fundamental parameter space of the models at GUT scale. It is then interesting to look back at 
the prospect of discovering gamma fluxes as a function of a physical parameter, the neutralino mass 
{iTi^), and to compare both SUSY breaking scenarios with mSUGRA. For that purpose, we have 
computed the total gamma flux as a function of m^ for three values of tan/3, scanning on the two 
other main parameters of each model ( [^^3/2, Sgs] in the moduli dominated case and [m^/2, 6+] in 
the dilaton dominated one). We have applied the experimental cuts coming from LEPII constraints 
on SUSY spectrum, higher order processes and neutralino relic density favored regions (the details 
are developed in appendix C). 

We have first plotted in Figs IBK and IHK the total gamma flux coming from the Galactic Center 
above 1 GeV as a function of the neutralino mass for a NFW profile and different values of tan /3 
(5, 20, 35) in the moduli and dilaton dominated scenarios. We note that: 
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FIG. 7: Dilaton-dominated scenario, a) The integrated gamma flux {E-, > 1 GeV) for a NFW profile in the dilaton 
parameter space as a function of the beta function 6+, for t — 0.25, p — 0, tan /3 = 5 (blue crosses) and 35 (red 
boxes) and b) the neutralino higgsino fraction as a function of b+ for the same values of tan 13. See the text for the 
comments. 



• For a fixed value of tan /3 and a given profile density, the regions covered by the two models 
do not overlap. In other words, if we obtain the value of tan/3 by another experiment, the 
measurement of the gamma-ray flux will allow to distinguish or even exclude SUSY breaking 
scenarios. 

• The flux is a decreasing function of the LSP mass, independently of the specific model 
adopted (see eq. I14|) : a heavier neutralino has a lower number density in the halo, and so a 
lower flux. Except some possible threshold effects, this is a general remark for any nature of 
indirect fluxes coming from neutralino annihilation at the Galactic Center. On the contrary, 
for lighter neutralino (dilaton-dominated scenario and mSugra) the threshold effect implies 
an increasing flux with rriy.. This effect is evaded for rriy. > mt, when the ti channel is open. 

• Fluxes in the dilaton breaking scenario are higher than in the moduli scenario (except for high 
tan/3 because of the j4-pole contribution, see below). The highest flux in moduli dominated 
case is obtained for low |(5gs| (Fig- EK)) where the neutralino is completely wino but does not 
have a sufficient relic density to fulfill the dark matter content of the universe (see Fig^. 
Those points are excluded by our cut on fi > 0.03. The remaining ones are situated in the 
near-zone "stau LSP" branch of Figs. HI and [T^ 

In the dilaton dominated breaking scenario, most of the points are excluded at low rriy. because 
they do not respect the accelerator constraints on m + (and m/j), whereas heavy neutralinos are 
forbidden by the closed parameter space limited by the gaugino condensation group /3-function 
coefficient 6+ (6+ < bEg ~ 0.57). This parameter space is larger for lower values of tan/3 because 
barely constrained by the realization of the electroweak symmetry breaking condition (^^ > 0). 
Most of the points with low integrated flux are excluded when imposing a higher limit on the relic 
density (fi^ < 0.3). This is due to the correlation between the fluxes and the relic density. In any 
case, if we take a lower upper omega limit (like the WMAP one) we keep the points with the best 
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FIG. 8: MODULI : The integrated flux of 7-ray above 1 GeV in a) and 60 GeV in b) per centimeter square per 
second versus neutralino mass for tan/? — 5, 20 and 35 taking a NFW profile for tlie lialo, after a scan in the moduli 
parameter space for t = 0.25, p = 0. The points that violate the accelerator constraints (see the appendix C for 
more details) are not represented here. Sensitivities of present and future experiments are represented in solid lines. 
We also show an "equivalent" mSUGRA cloud (see discussion in appendix B for the choice of parameter range). 
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FIG. 9: DILATON : The integrated flux of 7-ray above 1 GeV in a) and 60 GeV in b) per centimeter square 
per second versus neutralino mass for tan/3 — 5, 20 and 35 taking a NFW profile for the halo, after a scan in the 
dilaton parameter space for t — 0.25, p — 0. The points that violate the accelerator constraints (see the appendix 
for more details) are not represented here. Sensitivities of present experiments and projects are represented in solid 
lines. We also show an "equivalent" mSUGRA cloud (see discussion in appendix for the choice of parameter space). 



detection rates. As was pointed out in [l9|, the cosmologically favored region of the parameter 
space is similar to the Focus Point (FP) region present in mSUGRA [4]|, i.e. a mixed higgsino 
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FIG. 10: The 77 and 27 lines versus neutralino mass for tan/3 = 5, 20 and 35 taking a NFW profile for the halo, 
after a scan in the dilaton parameter space and moduli one for t — 0.25, p = 0, and a comparison with mSugra. The 
points that violate the accelerator constraints are not represented here (see the appendix for the cut applied). 



region leading to strong couplings and high gamma-ray fluxes dominated by ti {Z exchange) or 
WW (ZZ) final state processes. This explains the correlation between the dilaton and mSUGRA 
clouds on Fig. IHl 

In the moduli dominated scenario, the dependence on tan f3 is not only due to the mass spectrum 
or a larger parameter space, but to the annihilation processes governing the gamma-ray fluxes. 

Indeed, in the moduli breaking scenario, one can have a relatively light pseudo scalar A. The 

A — 
annihilation of the neutralino is thus dominated by the s channel XX ~^ ^b with a branching ratio 

approaching the 90 percents in most cases (the coupling Abb oc tan/?, gives a factor 50 in a between 

tan/J = 5 and 35). We can notice that thanks to the 6gs parameter the A-pole can be obtained 

for lower value of tan/3 than in mSUGRA. 

To study the observability of these fluxes, we have estimated the total gamma-ray flux above 
60 GeV, and compared it with the predicted sensitivity of the HESS experiment (Figs. |Hb andEb). 
The low neutralino mass region (< 400 GeV)is not accessible to EGRET but will be explored by 
GLAST, whereas HESS will give more precise information for a heavy neutralino. In that sense, 
GLAST and HESS experiments are complementary. Indeed, evading the NFW profile assumption 
and taking a less cuspy density, the GLAST accessible region of our models is much less overlapping 
the HESS one. For instance, assuming a NFW profile, a non discovery of any signal in HESS and 
GLAST will exclude any dilaton dominated scenario and the high tan/3 regime of the moduli 
domination. 

Finally, in the framework of effective heterotic models proposed in this paper, concerning direct 
and solar neutrino indirect detection, points satis fyin g all constraints in the dilaton case were all 
accessible to experiments whereas moduli was not [l9|. Concerning gamma indirect detection, the 
conclusion depends strongly on the halo profile assumption. For a NFW shape, the dilaton model 
is still the most attractive for detection but the moduli one can also be tested by experiment. In 
that sense, gamma indirect detection is less an explicit test of SUSY breaking scenario but a very 
complementary way to give some hints on particles physics hypothesis and dark halo astrophysical 
assumptions. 
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We close this section by showing 77 and Z7 hne results on Figs. 1101 Because dark mater 
in the halo is extremely non-relativistic, photons from these processes have an energy width of 
only AE'^/E'^ ~ 10~^ and are effectively mono-energetic. While this signal would be the most 
spectacular of all possible indirect signals, as we can see in Figs. El the processes being loop- 
suppressed are much too low to be detected, even in any string scenario studied here. On another 
hand, the integrated flux explored before is more observable, but far less distinctive and will 
certainly require additional confirmation to unambiguously distinguish it from the background or 
other exotic sources. 

V. SYNCHROTRON RADIATION 

Another interesting source for the indirect search of DM in the Galactic halo is the detection 
of the neutralino via the synchrotron radiation created by the propagation of the e^ products of 
neutralino annihilation, around the Galactic magnetic field. In our study, we have considered a 
magnetic field at equipartition at the Galactic center and constant elsewhere [38| 



B{r) = max 



5/4 
324^G ( — ) , 6^G 



(16) 



meaning that the field is constant for galacto-centric distances r > 0.23 pc. 

Lower values of the magnetic field imply a shift of the synchrotron spectrum toward lower 
energies. As a consequence, the flux at low frequency will increase, favoring the detectability of 
our models. Nevertheless, we prefer to be conservative and consider a B fleld at equipartition. 

The synchrotron flux per solid angle at a given frequency v is given by |2a] 

L,(V^)^^^f ^ "^'^" ^ ^,Y,{M,uy/^I{i^), (17) 

^^' 47r8V0.29^ e J M^ ^ ' ' ^^'' ^ ' 



where 



I(V)= / dsp2(r(s,V))i3-^/'(r(s,V)), 

J line of sight 



(18) 



and s is the coordinate running along the line of sight. Ye{M, v) is the average number of secondary 
electrons above the energy Em{i^) of the electron giving the maximum contribution at a given 
frequency z^ and given magnetic fleld B. We recall that 



4tt m^c^ u 



1/2 



In order to compare our predictions with the observational data, we to integrate over the 
corresponding solid angle, studying two distinct astrophysical situations and integrating over the 
string parameters of the theory in every case : 

• Flux at i^ = 408 MHz in a cone of half width 4 arc sec pointing around the Galactic Center 
in a NFW halo model. The observed flux is ~ 0.05 Jy 52]. 

• Flux aX, V = 327 MHz in a cone of half width 13.5 arc sec pointing around the Galactic 
Center in a NFW halo model. The observed flux is ~ 362 Jy |53| . 
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FIG. 11: Predicted synchrotron radiation from the Galactic center versus the neutralino mass for a NFW 
profile in a moduli dominated scenario and a dilaton one for tan/3 — 5, 20, 35 (after experimental cuts) at a frequency 
of 408 MHz (left) and 327 MHz (right) compared to the flux observed limit at the same frequencies (horizontal lines). 
We also show the results of a "corresponding" mSugra scan 

In our study, we have neglected two processes: self-absorption and absorption of electrons by 
the interstellar medium. In fact, it has been shown in Ref. |38| | that in a NFW profile, the optical 
depth r can be safely neglected unless very low frequencies are considered (of order ~ 1 Hz), 
which is not our case. Concerning the electron absorption, the authors of 'Sq observed that the 
absorption coefficient per length is such that a^ < 6 10~^^ pc~^{B/iJ.G){v/GHz)^^, justifying our 
approximation . 

In any case, for a NFW halo profile, the fluxes obtained for the models considered here do 
not constrain the mass of the DM particle, contrary to what happens in the case of Kaluza-Klein 
particles 38]. Naturally, the synchrotron flux shapes we obtain are similar to what we obtain for 
7 fluxes due to the same av dependence, but the resulting values in the case of a NFW profile 
are several orders of magnitudes lower than experimental constraints. Recently, the calculation 
of synchrotron radiation from neutralino annihilation has been revisited ||^]. Unfortunately the 
processes described there are not sufficient to make the flux in our models observable. 



VI. CONCLUSIONS 



We have studied gamma-ray and synchrotron emission from the Galactic center, in the context 
of an effective string inspired framework, and discussed two scenarios: in the first one, the SUSY 
breaking is transmitted by the compactification moduli T"; while in second one, it is transmitted by 
the dilaton field S, via their respective auxiliary fields. Typically, models in the dilaton dominated 
SUSY breaking scenario lead to a higher annihilation rate than the moduli scenario. Concerning 
the continuum gamma-ray flux, both scenarios are within the reach of the experimental sensitivities 
of GLAST and HESS for a NFW halo profile. For the same profile, the gamma-ray line signal is 
suppressed and beyond the experimental sensitivity. 

The synchrotron emission is too low to be constrained by experiments even with a more cuspy 
profile. Due to the dependence of the prospects of detection on both theoretical high energy physics 
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assumptions and astrophysical parameters of the halo dark matter distribution, gamma-ray indi- 
rect detection of neutrahno dark matter can give interesting information either on astrophysical 
hypotheses or on the SUSY breaking scenario by the complementarity with other experimental 
searches like direct detection or neutrino telescopes (studied in [1^ in the same context) or accel- 
erators. 
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APPENDIX A: ANNIHILATION BRANCHING RATIOS FOR tan/3 = 5. 

We show here the figures correspondmg to Fig0] and El for tan/3 = 5. The case of moduh 
domination regime is illustrated in Fig ll21 while in Fig. El we show the case of dilaton dominated 
SUSY breaking. 

Moduli : Branching ratios : tb=5 




FIG. 12: Moduli domination regime: dominant annihilation branching ratios in the ((5gs i "13/2) plane for 
tan/? = 5, t = 0.25, p = 0. Regions with gluino and stau LSP are indicated. We also show in grey the 
kinematically forbidden region for each channel. 
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Dilaton : Branching ratios : tb=5 




FIG. 13: Dominant annihilation branching ratios in the (6+,TO3/2) plane for tan/3 = 5, t — 0.25, p — 0. 
Regions where radiative electroweak symmetry breaking can not occur are indicated (No REWSB). We also 
show in grey the kinematically forbidden region for each channel. 



APPENDIX B: PARAMETER RANGES 



We performed a scan in the moduli and dilaton parameter space with the following values 

• Moduli : 

t = 0.25,p = 

tan/3 = 5 ; < my2 < 10000 GeV ; -30 < 5gs < 
tan/3 = 20 ; < mg/s < 16000 GeV ; -20 < 6gs < 
tan/3 = 35 ; < 7713/2 < 20000 GeV ; -10 < 6gs < 

• Dilaton: 

t = 0.25,p = 

tan/3 = 5 ; < mg/a < 14000 GeV ; < 6+ < 0.57 
tan/3 = 20 ; < 7713/2 < 4500 GeV ; < 6+ < 0.57 
tan/3 = 35 ; < 7713/2 < 4000 GeV ; < 6+ < 0.57, 
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FIG. 14: Low energy {m^,mi) and (Mi,/i) planes showing that the parameter space we choose are quite equivalent. 



corresponding to our previous study on direct detection and neutrino indirect detection [ig|. We 
also show for comparison an equivalent mSugra cloud with the following parameter : 

• mSugra : 
^0 = 

tan/3 = 5 ; < mo < 14000 GeV ; < m^/^ < 2000 GeV 
tan/3 = 20 ; < mo < 4000 GeV ; < mi/2 < 1500 GeV 
tan/3 = 35 ; < mo < 3000 GeV ; < mi/2 < 1000 GeV 

We show on Fig. 1141 two typical low energy planes: {my.,m^) and (Mi, /i) for tan /3 = 5, 35 only. 
This illustrates that our choices of parameter ranges are coherent from the low energy point of view 
and allow for comparison between the different models. An extended study of such comparison is 
the subject of a forthcoming paper J42l |. 



APPENDIX C: EXPERIMENTAL CONSTRAINTS 



We apply the following conservative cuts on our models: 
• Higgs mass: nih > 113.5 GeV |43l |. 



Chargino mass: m-^+ > 103.5 GeV J44| 



• Relic density: 0.03 < Q^h'' < 0.3, but we also show the WMAP [2^1 range O^j^^^^h 



0.1126 



-0.0161 
-0.0181' 



',WMAPr,2 



• b —>■ s^ Constraint |45|] : 

2.33 X 10"^ < BR(6 -^ sj) < 4.15 x 10-^ 

• The muon anomalous magnetic moment J4g]: 
JSM < 30.4 [2 a]. 



-11.6 



< 



rnew physics X'^^P 



23 



APPENDIX D: TOOLS 

For our computations we have interfaced SUSPECT !47(, MICROMEGAslii] and 
DARKSUSYlii. 

Concerning MSSM renormalization group equations and radiative electroweak symmetry break- 
ing we use the code SUSPECT. It includes 2-loops RGE evolution from high scale down to low 
energy scale (we took Q = \/i^2 tSfl]) and minimizes the 1-loop scalar potential by solving itera- 
tively the condition 



tan2 /? - 1 2' 



il' = ^^ °1 ,' : — ^Ml (Dl) 



giving the /i parameter which drives the crucial neutralino higgsino fraction for dark matter studies. 
We then transfer all MSSM parameters to the code MICROMEGAS which calculates the full 
physical mass spectrum including radiative corrections on Higgs masses and widths. The calcula- 
tion of neutralino relic density is then achieved following the iterative procedure described in [5l| 
and including all annihilation and co-annihilation processes by solving 

^ = y/ng,{T)/A5G{av) {Y^ - Y^q) (D2) 

where Y is the abundance at the temperature T, g^ a degrees of freedom coming from thermo- 
dynamics, G the Newton constant and {av) the thermally average cross section of all processes 
concerning Y . 

Masses and couplings are then entered in DARKSUSY which calculates the indirect detection 
rates of species i = 7,e^... coming from the Galactic Centre for a chosen galactic halo profile by 
splitting the particle physics dependance (oc (av) d(f) / dEi , (av) including here only annihilations) 
and the astrophysics part i.e the integration along the line of sight (J dsp'^{r{s, ip)), p{r) being the 
dark matter density distribution). 
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